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Abstract 

We add the Wess-Zumino-Witten term to the = 3 massive nonhnear sigma 
model and study the leading logarithms in the anomalous sector. We obtain the 
leading logarithms to six loops for vr*^ — )• 7*7* and to five loops for 7*7r7r7r. In 
addition we extend the earlier work on the mass and decay constant to six loops and 
the vector form factor to five loops. We present numerical result for the anomalous 
processes and the vector form factor. In all cases the series are found to converge 
rapidly. 
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1 Introduction 

Obtaining exact results in quantum field theory is rather difficult. One of the few things 
which can be easily calculated to all orders in renormalizable theories are the leading 
logarithms of the type {g"^ log /i^)" where fi is the subtraction scale and g the coupling 
constant. The analogue of this in effective field theories is not so simple since at each 
order in the expansion new terms in the Lagrangian appear and the recursive argument 
embedded in the renormalization group equations for renormalizable theories no longer 
applies. Nonetheless, one can calculate the leading logarithms in effective field theories 
using only one-loop calculations. This was suggested at two-loop order by Weinberg [1] 
and proven to all orders in [2]. 

In the massless case this has been used to very high orders in [3-5] for meson-meson 
scattering and form factors. In the massive case, many more terms contribute but in [6, 7] 
a method was developed for handling those, and the leading logarithms in the massive 
nonlinear sigma model were obtained to five-loop order for the mass, decay constant and 
the vacuum expectation value and to four-loop order for the vector and scalar form factors 
and the meson-meson scattering amplitude. A natural continuation of that program is 
to extend it to other sectors as well. We therefore add to the massive nonlinear sigma 
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model for = 3 the anomalous part via the Wess-Zumino-Witten (WZW) term. This 
allows us to study the leading logarithms for anomalous processes in two-flavour Chiral 
Perturbation Theory (ChPT). In addition one can hope that in this sector with its many 
nonrenormalization theorems it might be easier to guess the all order results when the first 
terms in the series are known. The WZW term only makes sense for = 3 so we do not 
work out the results for general N in this case. 

We have also improved the programs used in [6, 7] so that they now can be used to 
arbitrarily high orders given enough computing power. So, at least in principle, the problem 
of the leading logarithms is solved. In practice we obtained one order more than in the 
earlier work. 

The main part of the paper is devoted to the calculations of the leading logarithms 
(LL) for the two main anomalous processes in the pion sector, the full 71^7*7* and 7*7r7r7r 
vertices. For the former we obtained the LL to six loops and for the latter to five. The 
results indicate that in all cases the chiral expansion converges fast but we did not find a 
simple all-order conjecture. 

The results agree with all known relevant earlier calculations. As an additional check 
we have used several different parametrizations of the fields. 

In Sect. 2.1 we introduce shortly the massive nonlinear sigma model and in Sect. 2.2 
the two-flavour Wess-Zumino-Witten term. Sect. 3 contains the discussion of LL in the 
nonanomalous sector where we present our new results and show some numerical results. 
Here we also explain briefly the principles of the calculation. More details on the method 
can be found in [6, 7]. Sects. 4 and 5 are the main part of this paper. The LL are 
calculated in Sect. 4. We did not find a simple all-order conjecture but the LL indicate 
for example that the nonfactorizable part in 71^7*7* with both photons off shell should 
be small. We present numerical results in Sect. 5. In all cases we find good convergence. 
Sect. 6 shortly recapitulates our results. The appendix contains a dispersive argument to 
clarify the discrepancy with [4] for the vector form factor. 

2 The model 

2.1 Massive nonlinear 0{N + 1)/0{N) sigma model 

The 0{N + 1)/0(A^) nonlinear sigma model, including external sources, is given by the 
Lagrangian 

Cn. = ^D^^'^D^^ + F\^^ . (1) 

$ is a real A^-|-l vector, = $^ . . . $^), which satisfies the constraint = 1 and 
transforms under the fundamental representation of 0(A^ -|- 1). The covariant derivative is 

= 9^$" + t;"^ - . (2) 
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The vector sources are antisymmetric, = —vj^, and correspond to the unbroken group 
generators. The axial sources correspond to the broken generators. Lower-case Latin 
indices a,b, . . . run over 1, . . . , in the remainder and are referred to as flavour indices. 
The mass term contains the scalar, s^, and pseudoscalar, p", external sources as well 
as the explicit symmetry breaking term M^: 

= ( {2Bs' + ... p^) . (3) 

The vacuum condensate 

($^) = (10 ... 0) (4) 

breaks 0{N + 1) spontaneously to 0{N). We thus have in principle Goldstone bosons 
represented by 0. The explicit symmetry breaking term, the part containing M^, breaks the 
0{N + 1) symmetry to 0{N), enforcing the vacuum condensate to be in the direction (4) 
and gives a mass to the Goldstone bosons which at tree level is exactly M. 

This particular model is the same as lowest-order two-flavour ChPT for = 3 [8, 9] 
and has been used as a model for strongly interacting Higgs sectors in several scenarios; 
see, e.g., [10, 11]. 

The terminology for the external sources or fields is taken from two-flavour ChPT. The 
vector currents for N = 3 are given by v""^ = — ^^''^yc -^[^i^ ^abc ^.j^g Levi-Civita tensor. The 
electromagnetic current at lowest order is associated with v^. 

We write $ in terms of a real A^-component vector 0, which transforms linearly under 
the unbroken part of the symmetry group 0{N). We have made use of five different 
parametrizations in order to check the validity of our results. They are 
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$1 is the parametrization used in [8], $2 a simple variation. $3 is such that the explicit 
symmetry breaking term in (1) only gives a mass term to the field but no vertices. $4 is 
the parametrization one ends up with if using the general prescription of [12]. Finally, $5 
has been used by Weinberg in, e.g., [1, 9]. The reason for using different parametrizations 
is that for each one of them, the contributions are distributed very differently between the 
diagrams and obtaining the same result thus provides a thorough check on our calculations. 
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2.2 Wess-Zumino-Witten Lagrangian 

For = 3, the massive nonlinear 0{N + 1)/0{N) sigma model corresponds to two-flavour 
chiral perturbation theory, which is an effective field theory for QCD. It is well known 
that the chiral axial current of the latter is anomalous [13-16], leading to the occurrence 
of processes such as vr" — )■ 77 or ir'j — )■ nn. The Lagrangian that we have introduced above 
does, however, not account for this. The necessary interaction terms are contained in the 
Wess-Zumino-Witten term [17, 18], which must be added to the effective Lagrangian. It 
is constructed such that it reproduces the anomalous Ward identities. Kaiser derived the 
WZW term for two- flavour ChPT [19], where it is considerably simpler than in the case of 
three flavours. His result can be re-expressed in terms of the field $ as 

+ (9^$°$" - <l>°9^$°)<9pt;° + ^e"''^$°$^;>e5pt;°| . (6) 

The interaction with the axial current coming from the WZW term has been omitted. 
Note that the normalization of the vector field differs from the one used in [19]. The 
Lagrangian depends on the Levi-Civita tensor e"''"^ in the S0{3) flavour indices. This is an 
object specific to = 3. There is no obvious simple generalization^ to different so for 
anomalous processes we restrict the calculation to = 3. 

The Lagrangian in (6) is of chiral order implying that anomalous processes are of 
the same order at leading order, while one-loop corrections are already 0{p^). This is 
immediately clear from the presence of the epsilon tensor with four Lorentz indices: each 
one of them must be combined with either a derivative or an external vector field, both of 
which are 0{p). 

The interaction of the pseudoscalars with the photon field Afj, can be obtained from the 
WZW Lagrangian by setting the vector current to 

< = I A , K = eA,6'^'' ■ (7) 

The Lagrangian of (1) has also a symmetry that QCD does not have [18]. The fields 
under this extra symmetry, called intrinsic parity, transform as 

$o^$o^ a^^-a;, s^s, p^-p. (8) 

The Lagrangian (1) and higher orders are even under this symmetry while (6) is odd. 

As in the even sector, we have used several of the parametrizations given in (5). Intrinsic 
parity for the different is such that it is always odd, 0" — >■ —0*^. The WZW term leads 
to interactions of one, two, or three photons with an odd number of pions. The anomaly 
also generates purely mesonic interactions among an odd number of five or more Goldstone 

^The generalization to different n for the SU{n) x SU{n)/SU{n) case is easy but that is a different 
case than the 0{N + \)/0{N) considered here. 
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bosons. However, for two flavours the purely mesonic odd intrinsic parity processes vanish 
to all orders. 

We use anomalous and odd intrinsic parity as synonyms and refer to the = 3 case 
here occasionally as the two-flavour case since it corresponds to two light quark flavours. 



3 Leading logarithms in the even sector 

In this section we recapitulate and extend some results of the even intrinsic parity sector 
of [6, 7]. We focus on those results that will be needed for the later calculations in the 
anomalous sector. 

The leading dependence on log fi at each order, with /i the subtraction scale, is what we 
call the leading logarithm. It can in principle always be obtained from one-loop calculations 
as was proven using /3-functions in [2] and in a simpler diagrammatic way in [6]. We will 
discuss some of those results in the sections on the explicit calculation of the mass and 
77r — 7- vrvr. 

In effective field theories there is a new Lagrangian at every order. The observation 
of [6] that the needed parts of those Lagrangians can be generated automatically from 
the one-loop diagrams allowed to perform the calculations. The actual calculations were 
performed by using FORM [20] extensively. 

We have extended the programs used in [6, 7] so that they can in principle run to an 
arbitrary order^. The only limit is set by computing time, which grows rapidly with the 
order. The necessary additions include routines that generate the required diagrams at a 
given order and calculate one-loop diagrams with an arbitrary number of propagators. In 
this way, we have verified some of the earlier results and obtained the coefficient of one 
more order for the mass, decay constant and vector form factor. 

In effective field theories writing the expansion in terms of lowest order or physical 
quantities can make quite a big difference in the rate of convergence. We therefore follow [7] 
in using two different expansions in terms of leading logarithms. A given observable Ophys 
can be written in different ways: 

Ophys = Oo (1 + aiL + 02^2 + ■■■), (9) 
Ophys = Oo (l + ciLphys + C2Lphys H ) , (10) 

where the chiral logarithms are defined either from the lowest-order parameters M and F 
as 

or from the physical decay constant F.^ and mass as 



^phys 



^In [6, 7], a different program was used for each diagram with a given number of vertices/propagators 
and the list of possible diagrams was constructed by hand. 
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3.1 Mass 



In this subsection we will present the first few coefficients Oj and q of the expansions in (9) 
and (10) with Ophys = and Oq = M^. In addition, we have also calculated the generic 
two-point function, which is needed for the wave-function renormalization. We extend the 
result from [6] by giving the expansion coefficients and Cj also at sixth order. 

Let us briefly recapitulate the strategy that is followed to obtain the expansion coeffi- 
cients. The starting point is the Lagrangian (1) which generates vertices with an arbitrary 
even number of pion legs. These lowest-order vertices are diagrammatically denoted by \a\ 
with the corresponding number of legs appended. The zero in this symbol refers to the 
order in the expansion. If we want to calculate the leading logarithm for the two-point 
function at one-loop level, we must evaluate the tadpole diagram with one insertion of the 
leading-order four-pion vertex. Its divergence must be canceled by a counter term from 
the next order. This can be depicted schematically as 



(13) 



Also at the next order, divergences must cancel, but the situation is somewhat more compli- 
cated. There are two one-loop diagrams from which, using the results of [2, 6], the leading 
divergence can be determined. This thus determines the relevant part of the second-order 
Lagrangian: 






(14) 



We already have all information ready in order to calculate the second diagram: the tree- 
level vertex follows directly from the lowest-order Lagrangian and the next-to-leading-order 
vertex has been determined in (13). The first diagram, however, contains the next-to- 
leading-order vertex with four pion legs, which we do not know yet. In order to obtain its 
divergence, we must calculate two more diagrams: 



(15) 



The algorithm continues to higher orders in exactly the same way. All the diagrams needed 
for the two-point function up to third order are shown in [6], Figs. 3-5. The total number 
of diagrams needed for the mass up to order n is 1, 5, 16, 45, 116, 303, . . .. 
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i 


Oj for = 3 


a-i for general N 


1 


-1/2 


1-1/2 N 


2 


17/8 


7/4 - 7/4 iV + 5/8 


3 


-103/24 


37/12 - 113/24 N + 15/4 iV^ _ 


4 


24367/1152 


839/144 - 1601/144 N + 695/48 A^^ - 135/16 A^^ + 231/128 N'^ 


5 


-8821/144 


33661/2400 - 1151407/43200 N + 197587/4320 iV^ - 12709/300 iV^ 
+6271/320 - 7/2 


6 


1922964667 
6220800 


158393809/3888000 - 182792131/2592000 N + 1046805817/7776000 N"^ 
-17241967/103680 + 70046633/576000 N'^ - 23775/512 
+7293/1024 



Table 1: The coefficients of the leading logarithm L* up to z = 6 for the physical meson 
mass. 



i 


a for = 3 


Cj for general A^ 


1 


-1/2 


1 - 1/2N 


2 


7/8 


-l/4 + 3/4Af- 1/8A^2 


3 


211/48 


-5/12 + 7/24Af + 5/8 Af2 - 1/16 A^^ 


4 


21547/1152 


347/144 - 587/144 A^ + 47/24 Ar2 + 25/48 A^^ - 5/128 A^^ 


5 


179341/2304 


-6073/1800 + 32351/2400 A^ - 59933/4320 A^2 + 224279/43200 A^^ 
+761/1920 AT'^ - 7/256 A^5 


6 


2086024177 
6220800 


-17467151/3888000 - 10487351/2592000 A^ + 68244763/1944000 A^^ 
-5630053/172800 A^^ + 18673489/1728000 A^"^ + 583/2560 A^^ 
-21/1024 A^*' 



Table 2: The coefficients q of the leading logarithm -Z^phys up to i = 6 for the physical 
meson mass. 



We did not find a simple formula to estimate the total number of diagrams needed. We 
did find a conjecture, verified up to 12th order, about the number of diagrams needed with 
only two external legs at each order: 

/■ n— 3 

. ^ .^'^-2 + 3x2 2 -1 fornodd 

# two-pomt diagrams = ■{ n , (Id) 

~n-2 _^ 22 — 1, for n even 

that is: 1, 2, 4, 7, 13, 23, 43, 79, 151, . . . diagrams with two external legs at order n. 

The coefficients and q in the expansion of the physical mass are listed up to sixth 
order in Tables 1 and 2. The sixth-order results are new. We can use these results to check 
the expansions and how fast they converge. We chose F = 0.090 GeV, = 0.0922 GeV 
and yU = 0.77 GeV for the plots presented here in Fig. 1. 
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0.1 0.15 
[GeV^] 



0.1 0.15 0.2 

[GeV^l 



Figure 1: The contribution of the leading logarithms to M^/M"^ order by order for F = 
0.090 GeV, = 0.0922 GeV, /i = 0.77 GeV and = 3. The left panel shows the 
expansion in L keeping F fixed, the right panel the expansion in Lphys keeping F^, fixed. 



3.2 Decay constant 

The decay constant is defined by 

{0\jl,mp)) = ^F.P.S"' ■ (17) 

We thus need to evaluate a matrix-element with one external axial field and one incoming 
meson. The diagrams needed for the wave function renormalization were already eval- 
uated in the calculation for the mass in the previous subsection. What remains is thus 
the evaluation of all relevant one-particle-irreducible (IPI) diagrams with an external a^. 
At one-loop order there is only a single diagram, but at higher orders, we must calcu- 
late 2, 4, 7, 13, 23, . . . diagrams, not counting the auxiliary diagrams required for the 
renormalization of higher-order vertices with more than two legs. We do not show these 
diagrams here, but they can be found up to third order in [7]. The total number of dia- 
grams with an axial current that needs to be calculated for the decay constant to order n 
is 1,5,18,56,169,511,... 

We give the coefficients for both leading logarithm series with Ophys = F^^ and Oq = F 
in Tables 3 and 4. The sixth order is again a new result. Note that once the expression of 
-Ftt as a function of F is known one may express the remaining observables as a function 
of the physical and F^. This has already been used to calculate the coefficients q in 
Tables 2 and 4 from the corresponding a^. 

We have plotted in Fig. 2 the expansion in terms of the unrenormalized quantities and 
in terms of the physical quantities. In both cases we get a good convergence but it is 
excellent for the expansion in physical quantities. 
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i 


Cj for = 3 


tti for general N 


1 


1 


-1/2 + l/2iV 


2 


-5/4 


-l/2 + 7/8iV-3/8iV2 


3 


83/24 


-7/24 + 21/16 iV - 73/48 iV2 _^ 1/2 


4 


-3013/288 


47/576 + 1345/864 iV - 14077/3456 iV^ + 625/192 iV^ - 105/128 iV^ 


5 


2060147/51840 


-23087/64800 + 459413/172800 iV - 189875/20736 iV^ 
+546941/43200 iV^ - 1169/160 iV^ + 3/2 iV^ 


6 


69228787 
466560 


-277079063/93312000 + 1680071029/186624000 N 
-686641633/31104000 iV^ + 813791909/20736000 iV^ 
-128643359/3456000 iV^ + 260399/15360 iV^ - 3003/1024 iV^ 



Table 3: The coefficients Oj of the leading logarithm L* for the decay constant in the 
case = 3 and in the generic case. 



i 


a for iV = 3 


Cj for general 


1 


1 


-1/2 + l/2iV 


2 


5/4 


l/2-7/8iV + 3/8 7V2 


3 


13/12 


-1/24 + 13/16A^ - 13/12 iV2 _^ 5/15 iyS 


4 


-577/288 


-913/576 + 2155/864 - 361/3456 iV^ - 69/64 iV^ _^ 35/128 A^^ 


5 


-14137/810 


535901/129600 - 2279287/172800 A^ + 273721/20736 A^^ 
-11559/3200 - 997/1280 A^"^ + 63/256 A^^ 


6 


37737751 
466560 


-112614143/93312000 + 3994826029/186624000 A^ 
-1520726023/31104000 A^2 _^ 276971363/6912000 A^^ 

-39882839/3456000 A^^ - 979/15360 A^^ + 231/1024 A^^ 



Table 4: The coefficients q of the leading logarithm -^^phys fo'^ the decay constant in the 
case = 3 and in the generic N case. 



3.3 Vector form factor 

Before we proceed to the discussion of anomalous processes, we turn to the last ingredient 
from the even intrinsic parity sector which will be used later: the vertex involving a single 
photon and an even number of pions. It is directly connected to the vector form factor, 
which is defined by 



The procedure to find the leading logarithm for this observable follows the lines of the 
one for the decay constant. For the wave function renormalization one may again use the 
results obtained in the mass calculation. We express here the results in terms of t = t/M^ 
and 



log 



(19) 
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0.1 0.15 
[GeV^] 



0.1 0.15 
[GeV^l 



Figure 2: The contribution of the leading logarithms to F-,^/F order by order for /x = 
0.77 GeV and iV = 3. The left panel shows the expansion in L with F = 0.090 GeV fixed, 
the right panel the expansion in Lphys with F^^ = 0.0922 MeV fixed. 



with a scale Ai^ that is some combination of t and M^. Again we have added one more 
order compared to the result in [7]. To fifth order we find 



Fvit) = 1 + Lm 1/6 1 



+ L 



+ L3 



_v( t (-11/12 + 5/12 iV) +? (5/36 - 1/24 A^) 
t (+1387/648 - 845/324 iV + 7/9 iV^) + P (-4007/6480 
+ 3521/6480 iV - 29/180 A^^) + P (+721/12960 - 47/1440 iV + 1/80 iV^) 

+ t (-44249/15552 + 222085/31104 iV - 55063/10368 iV^ + 127/96 iV^) 

+ P (+349403/155520 - 15139/4860 iV + 86719/51840 iV^ - 199/480 7V^) 
+ P (-85141/155520 + 885319/1555200 TV - 5303/19200 iV^ + 21/320 iV^) 

+ t^ (+4429/103680 - 57451/1555200 iV + 289/14400 iV^ - 1/240 iV^) 



^ t (-2278099/777600 - 2377637/466560 iV + 64763783/4665600 iV^ 

- 178063/19200 + 69/32 N^) + P (-62212433/11664000 
+ 27685279/2332800 iV - 376597697/38880000 iV^ 

+ 53519593/12960000 iV^ - 361/400 iV^) + P (74033879/30240000 

- 2247054421/544320000 iV + 32125153/11340000 iV^ 

- 13264877/12096000 iV^ + 1209/5600 iV^) + P (-299603257/816480000 
+ 213192107/408240000 iV - 98330371/272160000 iV^ 

+ 546331/3780000 A^^ - 233/8400 7V^) + P (3090331/163296000 
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- 36097349/1632960000 iV + 28441883/1632960000 iV^ 

- 15971/2268000 iV^ + 1/672 7V^) . (20) 

Note that Fv{0) = 1 as it should be. 

The vector form factor was also calculated in the massless case in [4]. In order to 
transform our result into this limit, we define 

^*-™logf-^). (21) 



167r2F2 ° V t 

Replacing A^^ _i. f ^^d then performing the limit — t- (which implies F.,^ ^ F), we 
get 

F°(t) = l + Kt/6 + K^{5/36 - N/24) 
+/s:4^ (721/12960 - 47/1440 N + 

+/s:i^ (4429/ 103680 - 57451/1555200 iV + 289/14400 A^^ - iVV240) 
(3090331/163296000 - 36097349/1632960000 
+ 28441883/1632960000 - 15971 /2268000 A^^ + A^V672) . (22) 

This differs from the result of [1]. The difference is discussed in the appendix using a 
dispersive approach as an alternative check, which agrees with [7] and (22). Up to the 
given order, the coefficient of the highest power in in (22) at each order is of the form 

fv= , \ , hr ' (23) 



(n + l)(n + 2) V2 

such that, assuming this representation of fy to be valid also at higher orders, we can write 

oo 

F° (t) = l + J2 fv^ti^ + 0{l/N)) . (24) 

n=l 

The summation can be performed explicitly and we obtain the next-to-large A^ result in 
the chiral limit in a closed form: 



A^ KtN^ 



- ^. , KtN 

l-(l + ^)log(l' 



(25) 



These large A^ formulas are also present in [1] up to a sign mistake. In that article, the 
large A^ has been explicitly calculated to all orders. 

We close this section with giving the expansion for the radius and curvature of the 
vector form factor defined by 

Fvit) = l + Ur^)vt + cvt^ + ... . (26) 
6 

The coefficients Cj for the expansion in physical quantities are given in Tables 5 and 6 in 
units of M^, again adding one order compared to [7]. The result up to two- loop order 
agrees with the LL extracted from the full two-loop calculation [21]. We do not present 
numerical results for the vector form factor since these are dominated in the physical case 
A^ = 3 by the large higher-order coefficient contributions, see, e.g., [8, 21]. 
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i 


a for N = 3 


Q for general N 


1 
2 
3 

4 
5 


1 

2 

853/108 
50513/1296 
120401/648 


1 

-ll/2 + 5/2iV 

1387/108 - 845/54 iV + 14/3 TV^ 

-44249/2592 + 222085/5184 iV - 55063/1728 A^^ _^ 127/16 iV^ 
-2278099/129600 - 2377637/77760 iV + 64763783/777600 TV^ 
- 178063/3200 iV^ + 207/16 iV^ 



Table 5: The coefficients q of the leading logarithm Lphys in the expansion of the radius 
(r^)y in the case = 3 and for general A^. 



i 


a for TV = 3 


Cj for general 


1 








2 


1/72 


5/36 - 1/24 A^ 


3 


-71/162 


-4007/6480 + 3521/6480 AT - 29/180 A^^ 


4 


-25169/7776 


349403/155520 - 15139/4860 AT + 86719/51840 A^^ _ 199/48O A^3 


5 


-1349303/72900 


-62212433/11664000 + 27685279/2332800 A^ 
-376597697/38880000 A^2 + 53519593/1296OOOO A^^ _36l/400Ar4 



Table 6: The coefficients q of the leading logarithm Lphys in the expansion of the curvature 
cy in the case A^ = 3 and for general A^. 

4 Leading logarithms in the anomalous sector 

4.1 7r7 — > 7T7T 

The process 7r°7 — )■ tt^tt^ is forbidden by C-symmetry and we will therefore concentrate on 
7r~7 — 7- 7r^7r°. The latter can be represented by the anomalous VAAA quadrangle diagram 
at quark level. We follow the notation introduced in [22] where the one-loop order was 
calculated. At tree-level, the amplitude for 7r~(pi)7(fc) — )■ vr~(p2)7r°(po) can be obtained 
from the Wess-Zumino-Witten Lagrangian (6): 

Ao = ^^,e,.a0''{k)pMpl (27) 
For higher orders we express the results in terms of physical variables only as 

A = tF^^s, t, u)e^,„^£^(A;)KP2Po . (28) 
with the Mandelstam variables 

s = {pi + k)\ t = {pi-p2)\ u = {pi-poy, s + t + u = 3M^ + e . (29) 
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Figure 3: The irreducible diagrams for the process 717 — )■ tttt up to two-loop level. The 
first two diagrams are the one-loop level. 

The function F^'^{s,t,u) for 77r — t- tttt is fully symmetric in s,t,u. We write it in terms of 
Ft, as 

F-'^s,t,u)=F^-f{s,t,u), Ff = ^^. (30) 

n 

If one expands /(s, t, u) as a polynomial in s, t, m one can use the relation in (29) to see 
how many new independent kinematical quantities can appear at each order. Up to fifth 
order in s, t, u there is only one at each order. At sixth order there are two. For the first 
five orders we choose as independent quantities^ 

A„, = s" + r + . (31) 

We also define fc^ = fc^/M^ and A„ = An/M^*^. In the end we write f{s,t,u) in terms of 

^2,M2,A2,...,A5. 

The main focus of this article is calculating the leading logarithms in the anomalous 
sector. The procedure follows very similar steps as in the even sector. From the even sector 
we will need the wave function renormalization and the expressions for the higher-order 
purely mesonic vertices as well as the decay constant and mass when using an expansion 
in terms of physical quantities. Vertices coupling an even number of pions to a single 
photon are not needed at this point, because the two- flavour anomaly does not contain 
interaction terms involving an odd number of pions and no photon. What remains to be 
calculated are the irreducible diagrams with three external pions and one external photon. 
The required diagrams up to two-loop order are depicted in Fig. 3. As in Section 3, a box 
with n inside, [n], means a vertex of order n. To reach the one-loop level, we must calculate 
the first two of these diagrams. Inspection of the remaining diagrams for the two-loop level 
shows that all vertices are already known except the one with five pions and one photon 
in the third diagram. The diagrams that are needed in order to obtain its divergence are 
shown in Fig. 4. Note that the vertex with three pions and one photon in the sixth diagram 
of Fig. 3 is fixed by the one-loop calculation, the first two diagrams in the same figure. 

To go to higher orders we rapidly need vertices with many more legs. We have generated 
all diagrams needed and calculated them up to fifth order. We agree with the logarithm 

•^The same arguments can be applied to any process fully symmetric in s, t, and u. A prominent 
example is 77 — >• 37r°, where s + t + u = 3M^ + . 
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Figure 4: The irreducible (auxiliary) diagrams needed for the vertex bn^ up to one-loop 
level. 



determined from the full one-loop result [22]. The results were obtained in several different 
parametrizations with consistent results. 

We have calculated the amplitude for 7r~7 — )■ 7r~7r° up to five-loop level. In the general 
case, i.e., for A;^, ^ 0, we obtain 

n) = 1 + L^l (3 + P) + - 3)(P + 33) 



6 V / ■'"'72 
L^-^ (9OA3 - 640A2 - 8157 + 2105P + 8IP + + Lj^- 



, , , , , , , , I532A4 

1296^ ) '"'155520 

+A3(88538 + 1890P) - A2(577760 + 12240P + 540^) - 2433375 + 1296190p 

1 



-57430r + 480A;'' + 185A;' 



1} 

^326592000 



A5(13252156) 

-A4(160744570 + 518350p) + A3(1465187530 + 39593272p + 247260fc^) 
-A2(6756522937 + 257781206^ + 11188776^ - 9160P) - 6498695163 

+12675091794P + 801259373fc^ + 4780240p + 2948600P - 1832^°^ 



(32) 



The symmetry in s, t and u is obvious in this expression. Note that at second order A2 
does not appear even though it could be present a priori. 

We can check whether we find a simpler expression in the massless limit and for an 
on-shell photon, /c^ = 0. In this case we need to express the result in terms of the logarithm 

log , (33) 



where A is some combination of s, t, and u: 



/^^°(s, t, n) = 1 + - AgLi - A4Li + 3313039 ^ 

^ ^ ' ' ^ 72 ^ ^ 19440 ^ 81648000 ^ ^ ^ ' 

It is clear from symmetry considerations that there should be no term linear in La. We 
do however not know why the quadratic term is also absent. 
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Figure 5: First type of the irreducible diagrams contributed to vr" — )■ 77 up to three- loop 
order 

4.2 TT^ 77 

This is the most important process in the odd-intrinsic parity sector of QCD, since it is 
the process in which the chiral anomaly was discovered. It remains the main experimental 
test thereof. 

We started our discussion of anomalous processes with 717 — vrvr because the leading 
logarithms for this process could be calculated from our results in the even intrinsic parity 
sector and just one type of anomalous vertex. 

We define the reduced amplitude F^^^^ for 7r° — )■ j{ki)j{k2) 

A = e,,^p eT{h)el''{h) k'^k^, F^,,{kl kj) . (35) 

Fj^^^lkf, k"^) is symmetric under the interchange of the two photons. 

The irreducible diagrams contributing to 7r° — )■ 77 consist of two different types of 
one-loop diagrams. The first type contains the diagrams where both photons are attached 
to the same vertex, while in the diagrams of the second type, the two photons connect to 
two different vertices, only one of which is anomalous. We show here the diagrams needed 
for the calculation of the leading logarithms up to third order. Those of the first type 
are depicted in Fig. 5. There is one diagram at one-loop order, there are two at two-loop 
order, and four at three- loop order. The diagrams of the second type are depicted in Fig. 6. 
This time, there is one diagram at one-loop order, there are three at two-loop, and eight 
at three-loop order. The figures do not contain the auxiliary diagrams that are needed in 
order to determine the higher-order vertices with more than one pion leg. We only mention 
that up to three-loop order, one needs to calculate 11 and 23 diagrams for the first and 
second type of topologies, respectively. 

The first type of diagrams is a consistent subset if we only keep the terms with two 
vector sources in the Lagrangian (6). That this leads to identical results for several different 
parametrizations provides a thorough check on the correctness of our programs for this class 
of diagrams separately. 

We write the result with kf = kf/M^ in the form 

F^,,ikl kl) = ^^F,{kl)F,{kl)F,,ikl kl)F . (36) 
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Figure 6: Second type of the irreducible diagrams contributed to tt^ — )■ 77 up to three- loop 
order 

The functions F^{kf) and F^^{kl, k^) are defined to be equal to one for k"^ = 0. F^-y{kl, k^) 
contains only those parts that cannot be absorbed in the and thus gives the part 
that cannot be obtained as a product of single photon form factors. Finally F gives the 
corrections for the on-shell decay vr^ — )■ 77. 

We have calculated all contributions needed for the LLs up to six-loop order and ob- 
tained 

F = 1 - 1/6 + 5/6 + 56147/7776 + 446502199/11664000 
+ 65694012997/367416000 , 
F^ik"^) = 1 + L;vi(l/6 P) + Ll^{5/2A P + 1/72 P) + (71/432 P + 1/24 P 

+ 1/1296 P) + L^(-24353/31104 P + 4873/10368 P - 2357/31104 P 
+ 145/31104 P) + (-548440741/81648000 P + 9793363/3024000 P 

- 32952389/54432000 P + 487493/13608000 P - 2069/10886400 P°) 
+ (-3519465627493/102876480000 P 

+ 3560724235307/205752960000 P - 1524042680197/411505920000 P 
+ 4741599089/11757312000 P - 510932327/13716864000 P° 
+ 1775869/914457600 p2), 

F^^(p, A;^) = 1 + LlJilkl^ + L^^^Ai^ [ - 203/7776 + 29/10368 (kl + A;^) 

+ 1/216 (k{ + - 1/144 klkl] + L^fc^^a [ - 5983633/10206000 
+ 46103/1632960 {kl + kl) + 372113/11664000 (k{ + A;^) 

- 211/5443200 (kl + Ai^) - 394157/9072000 klkl - 4/25515 A;2^2(a;2 + kl)^ 
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+ L^kfkl [ - 1072421939773/205752960000 

+ 1444445383/6531840000 {kf + kl) + 10840553807/102876480000 (k^ + kt) 
+ 282016297/205752960000 (kl + kl) + 6157391/4115059200 (kl + kl) 

- 3852620057/29393280000 klkl - 154739/58320000 + kl) 

- 75041473/20575296000 klkl{kl + kl) + 174329/35721000 k\kl] . (37) 

The absence of the hnear term in F agrees with the statement from [23, 24] that the 
contribution from one-loop diagrams at NLO can be absorbed into F^. The quadratic term 
also coincides with the two- loop calculation of [25] and the complete one- loop expression 
for off-shell photons is the same as in [21]. 

Note that the nonfactorizable contribution F^^ only starts at three-loop order and that 
the part surviving in the chiral limit only starts at four-loop level. The leading logarithms 
thus predict this part to be fairly small. 

The lowest-order results for the two anomalous processes are connected via the current 
algebra relation [26, 27] 

F3-(0,0,0) = ^F^-,-,(0,0), (38) 

which holds exactly in the chiral limit. Even beyond this limit it is valid also at the one- 
loop level for the leading logarithms as can be seen by comparing (32) and (37). The 
current algebra relation remains true, if in both amplitudes one of the photons is allowed 
to be off-shell, i.e. 

F3-(s,t,u) =J^F (A;2o). (39) 



It turns out that in the chiral limit, this relation holds for the leading logarithms up to 
two loops, as can again be checked from (32) and (37). 



5 Phenomenology of the anomalous sector 

5.1 7r7 — > TTTT 

The only direct measurement of the 7r7 — tttt vertex has been performed at the IHEP ac- 
celerator in Serpukhov [28] using 7r~7 — 7r~7r°, where the 7 comes from the electromagnetic 
field of a nucleus via the Primakoff effect. The relevant cross-section was measured with a 
pion beam of = 40 GeV and the photons' virtuality was in the region /c^ < 2 x 10~^ GeV^, 
which can be neglected at the present precision. The analysis is for values of s < lOM^ 
and thus within the region where ChPT is applicable. Assuming the function F^'^{s,t,u) 
of (30) to be approximately constant, F^'^{s,t,u) ~ F^"^, they find 

FeJ, = 12.9 ± 0.9 ± 0.5 GeV"^ . (40) 
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Another derivation used the data on vr e ir e ir^ [29] to determine Fg^'^ [30]. They 
obtained 

F^-^p = 9.9 ± 1.1 GeV"^ or 9.6 ± 1.1 GeV^^ , (41) 

depending on the way electromagnetic corrections are included. 

The value in (40) needs to be corrected for extrapolation to the point s = t = u = 0. 
The one-loop ChPT corrections were evaluated in [22] and a possibly large electromagnetic 
correction identified in [31]. The main conclusion is that the experimental values are in 
fairly good agreement with the theoretical result of (30) which gives F^'^ = 9.8 GeV""^ as 
discussed further below. 

The comparison with [2S] goes via their result a/Z"^ = 1.63 ±0.23 ±0.13 nb and [28, 31] 

rf^ I o A/I 2 

min ymax 



^ J AMI \ Q'min ?max J ^ ~ 



iU247r 



1 . .... f 4M2^'/2 



f desin'^e\F^''f{s,t,u)\^ . (42) 
Jo 



Inserting the charged pion mass, q^^^^ = 2 ■ 10~^ GeV^, g^j^ = ((s — M^)/(2_E')) , and 
f{s,t,u) = 1 in (42) leads to the value in (40). 

The various known higher-order corrections can now be included via f{s,t,u): 

f{s, t,u) = l + /o™ ± ± /r + • • • • (43) 

The dependence on s,t,u is tacitly assumed for all functions /«. The index i refers to 
the h order. The leading-order electromagnetic correction f™ was determined in [31] as 
f™ = —2e'^F^/t and higher-order electromagnetic corrections were found to be small. 
The next-to- leading-order correction coming from one- loop graphs is, in the isospin limit, 
given by [22] 



/i 



loop ^ M! M2 



6F2 



(47r)2 V " /i 



l±31og^ + + 



(44) 



with J(s) = (s — 4M;^) J(s), where J(s) is the standard subtracted two-point function 

167r2J(g2) = ^log^-^ + 2 , a = Vl - 4M2/g2 . (45) 
cr ± 1 

The logarithmic term in (44) agrees with our LL calculation. The full expression accounting 
for the pion mass difference can be found in [31]. 

The contribution from the NLO Lagrangian can be expressed in terms of the low-energy 
constants introduced in [32] as 

/r^^ = 1287r'M^ (cf ± c^'-) . (46) 
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LO 


fEM 

Jo 


/r^(LL) 


rloop 
Jl 


(HLS) 


fl'''^ (CQM) 


fl'^^ (SDE) 


77l37r 

-^0 


12.9 


12.3 


12.0 


11.9 


11.4 


10.1 


12.0 



Table 7: The extraction of the anomalous 787? factor (in GeV~^) from experiment using 
various estimates of the higher-order corrections. LO includes no higher-order corrections 
and thus coincides with (40). The next column contains the QED correction /™. Then 
come the values including in addition the leading logarithm and the complete correction 
from The last three columns also contain /J;''*^'^ using the model estimates in (47). 

In order to estimate the value of fl'^'^'^, several methods exist in the literature: hidden local 
symmetry (HLS) [22], phenomenology [33], the constituent quark model (CQM) [33, 34], 
Schwinger-Dyson equation (SDE) [35], or resonance saturation [36], to name a few. The 
spread in results can be seen from the estimates following from three of these methods 

= = 0.048 (HLS), = 0.19 (CQM), = -0.01 (SDE) . (47) 
2Mp 

The two-loop corrections, estimated using dispersive techniques [37], are found to be small. 

These corrections can now be incorporated in the calculation of Fq'^ from the cross 
section measured at Serpukhov using (42). Turning them on one by one, we find the 
results listed in Table 7. Comparison of the third and fourth column shows that at the 
one-loop order, the leading logarithm provides a good estimate for the size of the complete 
correction: it accounts for 60% of the shift. The uncertainty on the listed values has two 
main sources: the experimental uncertainty of about ±1 GeV~^ and the model dependence 
of the cf^^. From the spread of the estimates in Table 7, the latter also is about ±1 GeV~^. 
The total error, adding quadratically, is thus about ±1.5 GeV~^, such that the theoretical 
result Fq'" = 9.8 GeV~^ agrees reasonably with the final values at the one-loop level. 

Let us now return to the discussion of the expansion /^^ in (32). The one-loop LL 
shifts the result by —0.3 as already shown. Adding the LL up contributions up to five-loop 
order in (42) leads to 

p3^LL ^ (^^2.9 - 0.3 + 0.04 + 0.02 + 0.006 + 0.001 + . . .) GeV"^ . (48) 

Clearly, the series converges rather well. The small size of the LLs beyond one loop indicates 
that the full corrections at higher orders are negligible. 

The total cross section obtained from only the LL contributions as a function of the 
center-of-mass energy is depicted in Fig. 7. 

5.2 7T^ — > 77 

For the decay tt^ — t- 77, there is more experimental information available. For a recent 
review, see [38]. The current PDG average ([39], updated 2011) for the lifetime of the 
neutral pion is based on six experiments: three relying on the Primakoff effect [40-42], a 
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Figure 7: The leading logarithm contribution to the total cross section for vr 7 — )> tt tt as 
a function of s. 

direct measurement [13], an e~^e~ collider measurement [11], and a measurement of the 
weak form factor in vr"'" — )■ e~^i"j [45], which is related to the tt'' lifetime via the conserved 
vector current hypothesis. This leads to the average lifetime r^o = (8.4 ± 0.4) x 10^^'^s. 
Including the recent precise measurement by PrimEx at JLab [46], 

r(7r° ^ 77)primEx = 7.82 ± 0.14 ± 0.17eV , (49) 

leads to a smaller uncertainty, t^^o = (8.35 ±0.31) x 10~^''s. Ongoing efforts by the PrimEx 
collaboration are expected to decrease the error by a factor two. 
The partial decay width is related to the decay amplitude by 

2 

■'^77 ~ 6471"^'^'^''"''^ ' ^^^^ 
At lowest order we find from the Wess-Zumino-Witten Lagrangian 

i^.^° = ^^ ^ r(7r°^77)L0^7.76eV, (51) 

which is in perfect agreement with the PDG average as well as with the PrimEx result (49). 
Higher-order corrections might still destroy the agreement and we can use the leading 
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logarithms to estimate the size of these contributions and to examine the convergence of 
the chiral series. Notice that no chiral logarithms are present at the one-loop level once 
everything is expressed in terms of the physical quantities F^^ and M^r [23, 24]. At the 
two-loop level, leading logarithms start to contribute [25]. At present the best prediction 
including electromagnetic and two- loop effects is [2-5] 



>77 



[8.09 ± 0.11) eV, (52) 



which leads to the lifetime r^o = (8.04 ± 0.11) lO^^'^s. 

Our result for F in (37) indicates that the convergence is fast and higher orders are 
small. Putting in /i = 0.77 GeV we obtain 

F = 1 + - 0.000372 + 0.000088 + 0.000036 + 0.000009 + 0.0000002 + . . . , (53) 

which clearly shows a fast convergence. 

We now turn to the discussion of the meson-photon transition form factor F^(— Q^), 
normalized to the value at = 0, which has been given in (37). It was measured by 
CELLO [47], CLEO [48], and recently by BaBar [49] mainly in the range 1 < < AO GeV^ 
New activity is expected for very low by KLOE-2 at DA$NE [50] which should directly 
test the prediction. 

The LL contribution up to fifth order has been given in (37). Our result for the LL 
contribution to this quantity is depicted in Fig. 8 together with the VMD prediction 



6 Conclusion 

In this paper we have extended the earlier work on leading logarithms in effective field 
theories to the anomalous sector. First we improved the programs used in the earlier work 
on the massive nonlinear sigma model [6, 7]. This allowed us to compute one order higher 
than in those papers and we presented results for the mass, decay constant and the vector 
form factor. For the latter we clarified the discrepancy with the chiral limit work of [4] 
and we presented some numerical results as well. 

The main part of this paper is the extension to the anomalous sector. We thus added the 
Wess-Zumino-Witten term to the massive nonlinear sigma model for = 3 and computed 
the leading logarithms to six-loop order for — )■ 7*7* and five-loop order for the 7*7r7r7r 
vertex. We did not find a simple guess for the coefficients which was one of the hopes when 
starting this work. In both cases the leading logarithms indicate that the chiral series 
converges fast and we presented some numerical results for the pion lifetime, the transition 
form factor and the 7*7r7r7r vertex. 
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Figure 8: The LL contribution to F^{—Q'^) at different orders. Also shown is the VMD 
prediction as a comparison. 
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A Dispersive approach for the pion form factor 

Since the leading logarithms for the vector form factor in the chiral limit obtained in [7] 
and here do not agree with the corresponding result from Kivel et al. [1], another check of 
our result is in order. In [5], it was found that the partial wave amplitudes for vrvr scattering 
are given by 

tl{s) = \Y.^L^^ In'^-^ (^) + (^(NLL), (55) 

n=l 
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with the dimensionless function S{s) = ^^^'^^ . The coefficients cj^; can be found in Tables I 
and II in [5]. 

The leading logarithms for the scalar form factor are 

Fs{s) = J2fnS{sr\n^'(-^] . (56) 

n=0 \ ^ / 

In this case, the results from [1] and [7] for the coefficients /(f are in agreement. The 
discontinuity across the cut of the scalar form factor must satisfy 

disc Fs{s) = t^.Fsis) , (57) 

which can be easily verified to hold for the coefficients given in [4, 7]. 
A similar expansion holds for the vector form factor: 

oo 2 

Ms) = Y.fnS{sr^n- ( - t) + ^(^LL) • 

n=0 

This time, however, the results from [ 1] disagree with ours and [7] for n > 2. The discon- 
tinuity across the cut of the vector form factor must hold 

disc Fv{s) = tlFvis) , (59) 

which is only given for the from us and [7]. We therefore conclude that this is the 
correct result. 

Dropping the factor (— 1)^"'"^ in (12) of [1] brings that result in agreement with ours. 
That there is indeed a misprint in [ 1] was conffimed to us by the authors and was stated 
in the PhD thesis of A. A. Vladimirov. 



References 

[1] S. Weinberg, Phenomenological Lagrangians, Physica A96 (1979) 327. Festschrift 
honoring Julian Schwinger on his 60th birthday. 

[2] M. Biichler and G. Colangelo, Renormalization group equations for effective field 
theories, Eur.Phys.J. C32 (2003) 427-442, [hep-ph/0309049]. 

[3] N. Kivel, M. Polyakov and A. Vladimirov, Chiral logarithms in the massless limit 
tamed, Phys. Rev. Lett. 101 (2008) 262001, [arXiv : 0809 . 3236]. 

[4] N. Kivel, M. Polyakov and A. Vladimirov, Leading chiral logarithms for pion form 
factors to arbitrary number of loops, JETP Lett. 89 (2009) 529-534, 
[arXiv: 0904. 3008]. 



23 



J. Koschinski, M. V. Polyakov and A. A. Vladimirov, Leading infrared logarithms 
from unitarity, analyticity and crossing, Phys.Rev. D82 (2010) 014014, 
[arXiv: 1004.2197]. 

J. Bijnens and L. Carloni, Leading logarithms in the massive 0(N) nonlinear sigma 
model, Nucl.Phys. B827 (2010) 237-255, [arXiv : 0909 . 5086]. 

J. Bijnens and L. Carloni, The massive 0(N) non-linear sigma model at high orders, 
Nucl.Phys. B843 (2011) 55-83, [arXiv: 1008.3499]. 

J. Gasser and H. Leutwyler, Chiral perturbation theory to one loop. Annals Phys. 
158 (1984) 142. 

S. Weinberg, Nonlinear realizations of chiral symmetry, Phys.Rev. 166 (1968) 
1568-1577. 

M. B. Einhorn, Speculations on a strongly interacting Higgs sector, Nucl.Phys. B246 
(1984) 75. 

F. Sannino, Dynamical stabilization of the Fermi scale: Phase diagram of strongly 
coupled theories for (minimal) walking technicolor and unparticles, 
arXiv: 0804. 0182. 

S. R. Coleman, J. Wess and B. Zumino, Structure of phenomenological Lagrangians. 
1., Phys.Rev. 177 (1969) 2239-2247. 

S. L. Adler, Axial vector vertex in spinor electrodynamics, Phys.Rev. 177 (1969) 
2426-2438. 

S. L. Adler and W. A. Bardeen, Absence of higher order corrections in the 
anomalous axial vector divergence equation, Phys.Rev. 182 (1969) 1517-1536. 

W. A. Bardeen, Anomalous Ward identities in spinor field theories, Phys.Rev. 184 
(1969) 1848-1857. 

J. Bell and R. Jackiw, A PC AC puzzle: — )■ 77 in the sigma model, Nuovo Cim. 
A60 (1969) 47-61. 

J. Wess and B. Zumino, Consequences of anomalous Ward identities, Phys. Lett. B37 
(1971) 95. 

E. Witten, Clobal aspects of current algebra, Nucl.Phys. B223 (1983) 422-432. 

R. Kaiser, Anomalies and WZW term of two flavor QCD, Phys.Rev. D63 (2001) 
076010, [hep-ph/00 11377]. 

[20] J. Vermaseren, New features of FORM, math-ph/00 10025. 



24 



J. Bijnens, G. Colangelo and P. Talavera, The vector and scalar form- factors of the 
pion to two loops, JHEP 9805 (1998) 014, [hep-ph/9805389]. 

J. Bijnens, A. Bramon and F. Cornet, Three pseudoscalar photon interactions in 
chiral perturbation theory, Phys.Lett. B237 (1990) 488. 

J. Donoghue, B. R. Holstein and Y. Lin, Chiral loops in , iff — )■ 77 and rj-rj' 
mixing, Phys. Rev. Lett. 55 (1985) 2766-2769. 

J. Bijnens, A. Bramon and F. Cornet, Pseudoscalar decays into two photons in chiral 
perturbation theory, Phys. Rev. Lett. 61 (1988) 1453. 

K. Kampf and B. Moussallam, Chiral expansions of the 7r° lifetime, Phys. Rev. D79 
(2009) 076005, [arXiv: 0901 .4688]. 

S. L. Adler, B. W. Lee, S. Treiman and A. Zee, Low-energy theorem for 
7 + 7^7r + 7r + 7r, Phys. Rev. D4 (1971) 3497-3501. 

M. Terent'ev, Process n"^ — vt^'tt^ in Coulomb field and anomalous divergence of 
neutral axial vector current, Phys.Lett. B38 (1972) 419. 

Y. Antipov, V. Batarin, V. Bezzubov, N. Budanov, Y. Gorin et. al.. Investigation of 
the chiral anomaly ^ Stt in pion pair production by pion in the nuclear Coulomb 
field, Phys. Rev. D36 (1987) 21. 

S. Amendolia, M. Arik, B. Badelek, G. Batignani, G. Beck et. al.. First measurement 
of the reaction tt^c — vr^vr^e, Phys.Lett. B155 (1985) 457. 

L Ciller, A. Ocherashvili, T. Ebertshauser, M. Moinester and S. Scherer, A new 
determination of the '-fir — t- tttt anomalous amplitude via 7r~e~ — t- tt'c'it^ data, 
Eur. Phys. J. A25 (2005) 229-240, [hep-ph/0503207]. 

L. AmetUer, M. Knecht and P. Talavera, Electromagnetic corrections to 
77r± 7r°7r±, Phys. Rev. D64 (2001) 094009, [hep-ph/0107127]. 

J. Bijnens, L. Girlanda and P. Talavera, The anomalous chiral lagrangian of order 
Eur. Phys. J. C23 (2002) 539-544, [hep-ph/0110400]. 

O. Strandberg, Determination of the anomalous chiral coefficients of order p^, 
hep-ph/0302064. Masters Thesis (Advisor: Johan Bijnens). 

J. Bijnens, The anomalous sector of the QCD effective lagrangian, Nucl.Phys. B367 
(1991) 709-730. 

S.-Z. Jiang and Q. Wang, Computation of the coefficients for p^ order anomalous 
chiral Lagrangian, Phys. Rev. D81 (2010) 094037, [arXiv: 1001 . 0315]. 



25 



[36] K. Kampf and J. Novotny, Resonance saturation in the odd-intrinsic parity sector of 
low-energy QCD, Phys.Rev. D84 (2011) 014036, [arXiv: 1104.3137]. 

[37] T. Hannah, The anomalous process ■jn — ?> hit to two loops, Nucl.Phys. B593 (2001) 
577-595, [hep-ph/0102213]. 

[38] A. Bernstein and B. R. Holstein, Neutral pion lifetime measurements and the QCD 
chiral anomaly, arXiv: 1112.4809. 

[39] Particle Data Group Collaboration, K. Nakamura et. al.. Review of particle 
physics, J.Phys.GG37 (2010) 075021. 

[40] V. Kryshkin, A. Sterligov and Y. Usov, Measurement of lifetime of the 7r° meson, 
Zh.Eksp.Teor.Fzz. 57 (1969) 1917-1922. 

[41] G. Bellettini, C. Bemporad, P. Braccini, C. Bradaschia, L. Foa et. al., A new 

measurement of the 7r° lifetime through the Primakoff effect in nuclei, Nuovo Cim. 
A66 (1970) 243-252. 

[42] A. Browman, J. DeWire, B. Gittelman, K. Hanson, D. Larson et. al.. The decay 
width of the neutral n meson, Phys.Rev. Lett. 33 (1974) 1400. 

[43] H. Atherton, C. Bovet, P. Coet, R. Desalvo, N. Doble et. al.. Direct measurement of 
the lifetime of the neutral pion, Phys.Lett. B158 (1985) 81-84. 

[44] Crystal Ball Collaboration, D. Williams et. al.. Formation of the pseudoscalars vr", 
rj and f]' in the reaction 77 — )■ 77, Phys.Rev. D38 (1988) 1365. 

[45] M. Bychkov, D. Pocanic, B. VanDevender, V. Baranov, W. H. Bertl et. al.. New 
precise measurement of the pion weak form factors in 7r+ — )■ e~^i"y decay, 
Phys.Rev. Lett. 103 (2009) 051802, [arXiv: 0804. 1815]. 

[46] PrimEx Collaboration, I. Larin et. al., A new measurement of the vr^ radiative 
decay width, Phys.Rev. Lett. 106 (2011) 162303, [arXiv : 1009 . 1681]. 

[47] CELLO Collaboration, H. Behrend et. al., A measurement of the vr^, 77 and 77' 
electromagnetic form factors, Z.Phys. C49 (1991) 401-410. 

[48] CLEO Collaboration, J. Gronberg et. al.. Measurements of the meson-photon 
transition form factors of light pseudoscalar mesons at large momentum transfer, 
Phys.Rev. D57 (1998) 33-54, [hep-ex/9707031]. 

[49] BABAR Collaboration, B. Aubert et. al.. Measurement of the 77* — ?■ vr'' transition 
form factor, Phys.Rev. D80 (2009) 052002, [arXiv: 0905. 4778]. 

[50] G. Amelino-Camelia, F. Archilli, D. Babusci, D. Badoni, G. Bencivenni et. al.. 
Physics with the KLOE-2 experiment at the upgraded DAcpNE, Eur.Phys.J. C68 
(2010) 619-681, [arXiv: 1003.3868]. 



26 



